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Procede de transmission d'un signal multiporteuse module avec 
offset (BFDM/OM) 

Le domaine de I'invention est celui de la transmission de signaux 
5 numeriques, bases sur des modulations multiporteuses. Plus pr6cis6ment, 

I'invention concerne la transmission, et notamment la modulation et la 
demodulation, des signaux multiporteuses biorthogonaux (BFDM/OM, en 
anglais « Biorthogonal Frequency Division Multiplex/ Offset Modulation »). 

Depuis plusieurs annees, les modulations multiporteuses ont suscit£ un 
10 grand interet. Celui-ci se justifie, en particulier, dans le cas des communications 

avec les mobiles, ou leur efficacite a deja ete demontree pour la diffusion des 
signaux radio avec, tout d'abord, le systeme « Digital Audio Broadcasting » 
(DAB, en fran^ais : « Diffusion Audionumerique ») [1] (par souci de 
simplification et de lisibilitd, toutes les references citees dans la presente 
15 description ont 6t€ regroupees en annexe E) mais egalement en transmission haut 

d6bit sur lignes bifilaires tdlephoniques avec les systemes ADSL (Asymmetric 
Digital Subscriber Line) et VDSL (Very high bit rate Digital Subscriber Line) 
[2]. 

Dans les schemas de modulations multiporteuses usuels, un ensemble de 
20 frequences porteuses, choisi de maniere a satisfaire des conditions 

d'orthogonalite en temps et en frequence, est multiplexe. C'est le systeme dit 
« Orthogonally Frequency Division Multiplex » (OFDM, en 
frangais : « multiplex de frequences orthogonales »). 

II est possible d'associer a chacune des porteuses une modulation sans 
25 offset (SM : « Synchronous Modulation ») ou avec offset (OM : « Offset 

Modulation »). On aboutit alors respectivement aux systemes OFDM/SM et 
OFDM/OM. En particulier, V association a chacune des porteuses d'une 
modulation d'amplitude en quadrature, sans ou avec « offset », produit, 
respectivement, les modulations OFDM/QAM (QAM, en anglais : « Quadrature 
30 Amplitude Modulation ») et OFDM/OQAM (OQAM : « Offset QAM »). Cette 
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derniere modulation fonctionne sans intervalle de garde et offre egalement une 
possibilite de choix plus etendue en ce qui concerne la fonction prototype [3], 
[4]. 

Toutefois, l'orthogonalite de l'OFDM ne lui assure 1'optimalite que dans 
5 le cas de canaux de transmission que Ton peut assimiler a un bruit additif blanc 

et gaussien. Dans tous les autres cas, 1'optimalite de l'OFDM n'est pas garantie. 

De ce point de vue, les modulations multiporteuses biorthogonales 
(BFDM) offrent des possibilites supplementaires et, en particulier, elles peuvent 
constituer un meilleur compromis vis-a-vis de canaux de type radio-mobiles qui 
10 sont a la fois dispersifs en temps et en frequence [5]. 

Par ailleurs une modulation biorthogonale avec offset (BFDM/OM) 
permet de conserver l'avantage de l'OFDM/OM avec la possibilite d'obtenir des 
fonctions prototypes bien localisees en temps et en frequence. 

A titre indicatif, on rappelle brievement, en annexe A, les definitions 
15 essentielles concernant les aspects mathematiques lies aux modulations de type 

BFDM/OM. Ces aspects ont deja fait l'objet de publications, avec l'appelation 
BFDM/OFDM, egalement conservee dans les annexes de la presente description. 

Dans un article recemment sounds [6], une technique de discretisation des 
systemes de modulation BFDM/OM a deja ete proposee. Toutefois l'approche 
20 decrite en [6], [7] se base essentiellement sur la discretisation des equations 

continues qui etendent en discret le formalisme introduit en continu dans la 
reference [4] pour l'OFDM/OM. 

Pour rOFDM/OM, cela suppose done l'utilisation d'une transformation 
mathematique, puis de la transformation inverse (classiquement FFT 1 puis FFT). 
25 On tronque ensuite le signal discretise. 

L* invention a notamment pour objectif de fournir une nouvelle technique 
de modulation et de demodulation d'un signal BFDM/OM qui soit plus efficace 
et plus aisee a mettre en ceuvre que les techniques connues. 

Ainsi, un objectif de 1' invention est de fournir de telles techniques de 
30 modulation et de demodulation permettant d' assurer, sur un plan theorique, que 
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TIES (Interference entre Symboles) et 1'IEC (Interference entre Canaux) soient 
exactement nuls, sur un support fini. 

LMnvention a egalement pour objectif de fournir de telles techniques qui 
permettent de realiser des dispositifs satisfaisant structurellement Fannulation de 
5 VIES et de FEEC. 

Un autre objectif de Finvention est de fournir de telles techniques, 
permettant aussi bien la mise en ceuvre de fonctions prototypes sym&riques ou 
non, et identiques ou non a 1'emission et & la reception. 

Encore un autre objectif de Finvention est de fournir de telles techniques 
10 de modulation et de demodulation, qui permettent de reduire et de controler les 

retards de reconstruction, par exemple pour des applications temps reel ou 
interactive. En d'autres termes, un objectif est de fournir de telles techniques 
permettant, pour des filtres prototypes de longueur donnee, d'obtenir des ddlais 
de reconstruction qui ne sont pas fixes (et qui peuvent done etre plus faibles que 
15 ceux de FOFDM/OM). 

LMnvention a egalement pour objectif de fournir de telles techniques, qui 
soient optimales, par rapport k des distorsions produites par un canal gaussien 
et/ou par des canaux non gaussiens qui ne se r^duisent pas simplement a un bruit 
additif blanc et gaussien. 
20 Encore un autre objectif de Finvention est de fournir de telles techniques, 

permettant d'obtenir des performances superieures aux techniques connues, en 
terme de localisation de la transformee. 

L' invention a egalement pour objectif de fournir des dispositifs de 
modulation et/ou de demodulation, et plus generalement de transmission et/ou de 
25 reception de signaux, qui soient ais6s et peu couteux a realiser et a mettre en 

ceuvre. 

Ces objectifs, ainsi que d'autres qui apparaitront par la suite, sont atteints 
selon Finvention k l'aide d'un procede de transmission d'un signal multiporteuse 
biorthogonal BFDM/OM, mettant en oeuvre une structure de transmultiplexeur 
30 assurant : 
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une etape de modulation, a 1' aide d'un banc de filtres de synthese, 

presentant 2M branches paralleled M > 2, alimentees chacune par 

des donnees source, et comprenant chacune un expanseur d'ordre 

M et des moyens de filtrage ; 
5 - une etape de demodulation, a l'aide d'un banc de filtres d'analyse, 

presentant 2M branches parallels, comprenant chacune un 
decimateur d'ordre M et des moyens de filtrage, et delivrant des 
donnees recues representatives desdites donnees source, 
lesdits moyens de filtrage etant deduits d'une fonction de modulation prototype 

10 pr6determinee. 

En d'autres termes, l'invention propose une realisation nouvelle des 
systemes de modulation BFDM/OM, bas6e sur une description nouvelle d'un 
systeme de modulation, sous la forme d'un transmultiplexeur, appel6 par la suite 
transmultiplexeur module. Comme cela apparaitra par la suite, cette technique 
15 pr6sente de nombreux avantages, tant en termes de modes de realisation que 

d'efficacite" des traitements, et notamment de 1'annulation de 1TES et de 1'IEC. 

On notera qu'une telle structure de transmultiplexeur module, 
permettant la transmission d'un signal multiporteuse moduli avec offset, se 
distingue fortement des structures de transmultiplexeurs de l'art anterieur. En 
20 effet, les schemas connus de transmultiplexeurs presentent des facteurs de 

decimation-expansion inferieurs ou egaux au nombre de sous-bandes mises en 
oeuvre-L'approche de l'invention permet en revanche, par la mise en ceuvre, sur 
chacune des branches des bancs de filtres, de moyens de filtrage deduits d'une 
fonction de modulation prototype pred6terminee, d'obtenir un nombre de sous- 
25 bandes superieur (double) au facteur d' expansion et de decimation. 

De plus, une telle structure de transmultiplexeur module selon 
l'invention pnSsente l'avantage, par rapport aux transmultiplexeurs de l'art 
anterieur, de permettre un grand choix de filtres prototypes. 
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Preterentiellement, lesdits moyens de filtrage dudit banc de filtres de 
synthese el/ou dudit banc de filtres d' analyse sont respectivement regroupes 
sous la forme d'une matrice polyphase. 

Cela permet, sur un plan pratique, de simplifier la complexity 
5 operatoire du transmultiplexeur. 

De fagon avantageuse, au moins une desdites matrices polyphases 
comprend une transformee de Fourier inverse h 2M entrees et 2M sorties. Les 
inventeurs ont en effet montre que Tutilisation d'une telle transformee, pour 
laquelle des algorithmes sont disponibles (IFFT), permet de simplifier fortement 
10 la realisation et la mise en ceuvre de 1'invention. 

L'invention concerne egalement le procede de modulation d'un signal 
transmis selon le procede de transmission decrit ci-dessus. Un tel procede de 
modulation met avantageusement en ceuvre une transformee de Fourier inverse 
alimentee par 2M donnees source ayant chacune subie un decalage de phase 
15 predeterminee, et alimentant 2M modules de filtrage, suivis chacun d'un 

expanseur d'ordre M, dont les sorties sont regroupees puis transmises. 

L'algorithme de modulation peut alors delivrer des donnees s[k] telles 



20 



que 



2M-1 In k D-M In 

x}(n) = j2^x° k (n)e ™ 2 e 2M 

m-l 



= + 2kM)x\{n - 2k) 



k=0 



25 4*] = 



oix D = aM - (3, 

avec a entier presentant le retard de reconstruction ; 
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(3 entier compris entre 0 et M-l ; 
et 1_.J est la fonction "partie entiere". 

De la meme facon. Invention conceme le precede de demodulation 
d'un signal transmis selon le precede de transmission decrit precedemment. Ce 
precede de demodulation met avantageusement en ceuvre une transformee de 
Fourier inverse alimentee par 2M branches, elles-memes alimentees par ledit 
signal transmis, et comprenant chacune un decimateur d'ordre M suivi d'un 
module de filtrage, et alimentant 2M multipliers de decalage de phase, 
delivrant une estimation des donnees source. 

Le precede de demodulation peut ainsi, avantageusement, delivrer des 

donnees telles que : 

xf(n-a) = s[nM-P-l] 

Jc;'(n -«) = £>(/ + 2kM)x?(n -a -2k) 



k=0 



2tz D+M 2M-1 j—kl 
k=0 



. 2jt ,D+M 



=2M&"' , "' 2 IEPT(i,"(n-a), ".^«-,('"-«)X'l 

De facon avantageuse, dans le precede de modulation et/ou de 
demodulation lesdits modules de filtrage sont realises sous l'une des formes 
appartenant au groupe comprenant : 

les filtres a structure transverse ; 
les filtres a structure en echelle ; et 
les filtres a structure en treillis. 
D'autres structures de filtres peuvent bien sur etre envisagees, et 
notamment les structures de filtres a reponse impulsionnelle infinie (RID- 
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Selon un mode de realisation particulier, correspondant notamment k la 
structure en treillis, ledit signal multiporteuse biorthogonal est un signal 
OFDM/OM. Des solutions techniques particulieres peuvent alors etre 
envisagees. 

L'invention concerne egalement, bien sur, les dispositifs d'emission 
et/ou de reception d'un signal BFDM/OM, mettant en oeuvre les proc6d£s 
presentes ci-dessus. 

D'autres caracteristiques et avantages de T invention apparaitront plus 
clairement h la lecture de modes de realisation preferentiels, donnes a titre de 
simples exemples illustratifs et non limitatifs, et des dessins annexes parmi 
lesquels : 

la figure 1 illustre la structure generale d'un transmultiplexeur 
associe a la modulation BFDM/OM, selon Tinvention ; 
la figure 2 presente, de fa9on simplifee, une vue globale de la 
chaine mettant en ceuvre un transmultiplexeur tel qu' illustre en 
figure 1 ; 

la figure 3 est une representation sous une forme polyphase du 
transmultiplexeur de la figure 1 ; 

la figure 4 illustre sur un cas elementaire, V insertion d'un 
operateur de retard plac6 entre un expanseur et un decimateur, 
utilise dans la mise en ceuvre de la representation polyphase de la 
figure 3 ; 

les figures 5 et 6 presentent respectivement un modulateur et un 
demodulateur BFDM/OM realises a Taide d'une FFT inverse ; 
les figures 7 et 8 presentent des filtres a structure en echelle 
pouvant etre utilises a la place des filtres polyphases des figures 5 
et 6, respectivement lorsque s, parametre entier defini par la suite, 
est pair ou impair ; 
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la figure 9 illustre une structure sous la forme de treillis pour les 
filtres polyphases des figures 5 et 6, dans le cas d'un signal 
OFDM/OM a filtres prototypes symetriques ; 
la figure 10 presente un treillis selon la figure 9, dans le cas 
normalise : 

les figures 11A et 11B d'une part et 12A et 12B d'autre part 
illustrent les reponses temporelle et frequentielle obtenues dans 
deux modes de realisations particuliers, correspondant aux 
tableaux de l'annexe D. 
Comme indique precedemment, la technique de Invention repose 
notamment sur une approche particuliere de discretisation, visant a obtenir 
directement une description du systeme de type transmultiplexeur module. Outre 
l'avantage d'un cadre de description plus general, cette approche offre de 
nombreuses possibility d' exploitation des liens entre les bancs de filtres et les 
transmultiplexeurs, pour l'optimisation des structures de realisation et du calcul 

des coefficients associes. 

Apres avoir pr6sente la structure generate de representation des systemes 

BFDM/OM sous la forme d'un modele discret de type transmultiplexeur, on 

presente ci-apres quatre modes de realisation particuliers de l'invention 

correspondant respectivement a : 

deux modes de realisation BFDM/OM qui, au modulateur et au 
demodulateur, utilisent tous deux un algorithme rapide de 
transformee de Fourier inverse (IFFT) et se distinguent par le type 
d'implantation des composantes polyphases du filtre prototype : 

- Mode 1 : Algorithme IFFT + filtrage polyphase transverse ; 

- Mode 2 : Algorithme IFFT + filtrage en echelle. 

- deux modes de realisation adaptes a l'OFDM/OM, deduits du 
BFDM/OM : 

- Mode 3 : variante du Mode 1 verifiant l'orthogonalite discrete 
de l'OFDM/OM avec filtrage polyphase transverse et 
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possibility de mise en oeuvre de filtre prototype symetrique ou 
non ; 

Mode 4 : variante du mode 2 verifiant P orthogonality discrete 
de TOFDM/OM avec filtrage polyphase realise par une 
structure en treillis. 
Des mdthodes de design de filtres prototypes illustrant ces proc6des de 

realisation des modulations BFDM/OM et OFDM/OM sont dgalement 

presentees. 

Les resultats presentes illustrent notamment : 

les possibilitds supplementaires du BFDM/OM pour lequel le 
d61ai de transmission reste modulable pour une longueur de filtre 
prototype donnee. Ceci permet, par exemple a retard de 
transmission identique, d'ameliorer les performances en terme de 
localisation temps -frequence de la transformation associee au 
modulateur. Cela permet egalement de maintenir des 
performances elevees du point de vue de la selectivity tout en 
reduisant le delai de transmission ; 

dans le cas de systemes dits dos-a-dos, la possibility avec les 
modes 2 et 4 d'annuler totalement 1* interference entre symboles 
(IES) et F interference entre canaux (IEC) et d'obtenir ainsi ce que 
Ton peut egalement appeler la reconstruction parfaite. 
D f autres exemples non reportes ici, montrent Egalement qu'il est possible 
d'obtenir, en biorthogonal, des performances en localisation comparables k celles 
de TOFDM/OM, et ceci avec des filtres prototypes beaucoup plus courts. 

Pour faciliter la lecture, on retient les notations suivantes : les ensembles, 
par exemple R le corps reel, ainsi que les vecteurs et matrices, par exemple E(z) 
et R(z) les matrices polyphases, sont notes en caractere gras. Sinon Tensemble 
des symboles mathymatiques utilises est note en caractere standard avec, en 
general les fonctions du temps en minuscules et les fonctions des domaines 
transformys (z et Fourier) en majuscules. 
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1 - Formulation sous form * d'un transmnltiplexenr moduli 

A partir d'un filtre prototype causal p[k], deduit de h(t) par translation et 
discretisation, nous obtenons un schema de realisation qui est celui de la figure 1. 

Dans ce sch6ma les filtres F,{z)\\ et H,(z) 12, avec 0 < f < 2M - 1, se 
deduisent de p[k] (ou P(z)) par modulation complexe. aet fr0< P<M-1, sont 
deux entiers qui se relient a un parametre D de la modulation D = aM - /3. Les 
calculs permettant d'aboutir a ce schema sont reportes en annexe B. 

On peut noter egalement que les filtres prototypes peuvent etre different*. 
Par la suite on se contentera d'&udier le cas particulier ou q[k] = p[D - k], sans 
que cela ne limite la portee de la demande de brevet. 

La realisation d'un schema de modulation et de demodulation 
directement selon cette figure 1 serait extrSmement couteuse, en termes de 
complexite operatoire. Selon l'approche de l'invention, on decompose done les 
filtres prototypes P(z), en fonction de ses composantes polyphases G{z), ainsi 
que cela est presente en annexe C. 

Cette annexe C precise egalement la relation d'entree-sortie, les 
conditions a respecter sur les composantes polyphases et le retard de 
construction. 

7. - Rxemples de realisation 

L'ensemble des modes de realisation decrits par la suite est bas6 sur la 
mise en ceuvre d'une transform6e de Fourier discrete (TFD). 

Cette technique presente bien entendu 1'avantage que la TFD se traduit 
par des algorithmes de calcul rapides, intitules selon leur sigle anglo-saxon FFT, 
ou EFFT pour la transformee inverse. (On notera que les equations rdferencees 
(1) a (54) se trouvent dans les annexes A a C). 

On note : 

(\ o 



w; =%/2 



. 2ir D+M 
J 2M 2 



e 



(55) 
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(56) 



Et IV la matrice de la transformee de Fourier discrete de taille 2M x 2M : 



[W] kt = e 2M , 0<U<2M-1 



(57) 



En utilisant les equations (35) a (38) (annexe C) on obtient 



10 



( 0 



{G 0 (z 2 ) 



0 



E{z 2 ) = W[W' 



W*W 2 



(58) 



r G 0 (z 2 ) 0 "1 



(59) 
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On en deduit les schemas du modulateur de la figure 5 et du 
demodulateur de la figure 6, tous deux realises a l'aide d'une transformee de 
Fourier inverse IFFT 51, 61. Sur ces figures 5 et 6, s est un entier d6fini par 
D = 2.S.M + d, d 6tant un entier compris entre 0 et 2M-1, 

Les notations et donnees apparaissant sur les figures 5 et 6, de raeme que 
sur les autres figures, font bien sur partie integrante de la presente description. 

Pour simplifier, mais sans perte de generality, on suppose par la suite que 
le filtre prototype P(z) est de longueur 2mM de sorte que toutes les composantes 
polyphases sont de meme longueur m. 

2.1 - Mode 1 .Algo rithme IFFT et decomposition polyphase 

En reprenant les notations des figures 5 et 6, on deduit les algorithmes de 
modulation et de demodulation suivants, deja mentionnes plus haut : 
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2.7.7 Algorithme de modulation 

-v» = a m y ? " (60) 



J—- k — W (61) 



= 2MV2IFFl|x 0 0 (n),---,A- 2 V,(»)e y2 *'" 1 JW (62) 



•xf («) = Xp(' + 2 Jkiltf) Jti (n - 2 A) 



14] 

s[k] = 



(64) 



10 2.1.2 Algorithme de demodulation 

x']in-a) = s[nM-f3-l] ( 65) 

j?J (« - a) = §p(/ + 2AM)Jc , / 2 (n-a- 2*) (66) 

. In: D+M 2M-\ j—kl 



.2* ,D+M 



= 2M-V2V 2W ' 2 IFFT^," (n - a). - *ui-i<» ~ «))M 



20 « m . B -a=^ 



-f-^.o ( „_ a) | 



(68) 



(69) 



2.2- Mode 2 : IFFTet structure en echelle 

Les schemas en echelle constituent un moyen d' implantation propose 
recemment pour la realisation des bancs de filtres. Les inventeurs ont valide 
25 mathematiquement leur application a la BFDM/OM, decrite ci-apres. 
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On a vu qu'on peut ecrire une modulation BFDM/OM sous la forme d'un 
transmultiplexeur utilisant deux FFT inverses (figures 5 et 6), dans lequel 
apparaissent explicitement les composantes polyphases du prototype utilise. 
Chaque filtre polyphase peut alors s' ecrire sous la forme d'une 6chelle. Selon 
que s est pair ou impair, on peut remplacer les filtres G,(z) des figures 5 et 6 par 
les schemas donnSs par les figures 7 et 8. 

Pour aboutir k de tels schemas, une decomposition matricielle des 
composantes polyphases est mise en ceuvre, qui s'appuie sur des matrices 2x2 
dont le nombre et la nature se determinent en fonction de la longueur du 
prototype et du delai de reconstruction souhaites. 

Par exemple, pour generer la paire de composantes polyphases 
[G/(z),G w+ jfe)] nous procedons en deux etapes : 



T initialisation est rdalisde par un couple (Fq, F,). F 0 correspond k 
un produit de trois matrices auquel vont correspondre les trois 
premiers elements des schemas du haut aux figures 7 et 8. La 
forme exacte de Fj depend de la parite du parametre s. C'est la 
matrice identite pour s pair (cf. schema du haut a la figure 7) ou 
c'est un produit de deux matrices C 0 et B 0 pour s impair, auquel 
vont correspondre les deux elements suivants dans le schema du 
haut a la figure 8. On obtient ainsi un prototype de longueur 2M 
(s = 0) ou4M(5=l). 

Pour augmenter la longueur de P(z), sans ou avec accroissement 
de delai, on applique ensuite un jeu de matrices qui sera 
respectivement soit (A, B\ soit (C, D). On obtient ainsi la suite du 
schema de realisation. 



Le meme principe s' applique aux composantes polyphases 
[G d _f(z) 9 en prenant cette fois les inverses des matrices precedentes. 

Un interet de cette structure est qu'elle garantit une reconstruction 
parfaite, meme en presence d'erreur sur les coefficients calculus, en particulier 
des erreurs de quantification. 
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Par ailleurs, cette structure facilite aussi 1' optimisation du filtre prototype, 
par exemple en prenant un critere de localisation ou de selectivity en frequence : 

-j + lj(2/n + l) = mM coefficients au lieu de 2mM, 



M-. 



il suffit d'optimiser y 

sans introduire de contrainte de reconstruction parfaite. 
% - Complevite. des diffe rent^ realisations 

Pour effectuer une comparaison des differents modes de realisation 
proposes, on se place dans le cas commun oixN = 2mM. Dans ce cas, chaque 
composante polyphase a une longueur egale a m. 

Chaque composante polyphase peut Stre realisee sous forme transverse, 
sous la forme d'une echelle ou, dans le cas orthogonal, sous celle d'un treillis. 
Meme si les echelles et les treillis possedent deux sorties, une seule est 
exploitable. 

Sur chaque sous-bande, on effectue au niveau du modulateur : 

- une pre-modulation (un decalage de phase, c'est-a-dire une 
multiplication complexe) ; 

- une transformed de Fourier inverse ; 

- un filtrage polyphase. 

Au demodulateur, on effectue les memes operations dans le sens inverse. 

On peut done en deduire la complexite du transmultiplexeur complet 
avec pr6-modulation, en terme d'operations complexes (tableau 1) ou reelles 
(tableau 2). 





Additions complexes 


Multiplications 
complexes 


Realisation transverse 
Realisation en echelles 
Realisation en treillis 
(normalise) 


2m-2 + 2log 2 2M 
4m + 2 + 2 log 2 2M 
4m^4 + 21og 2 2M 


2m + 2 + log 2 2M 
4m + 2 + log 2 2M 
4m + 2 + log 2 2M 



TABLEAU 1 - Nombre d'operations complexes par sous-bande et par 
echantillon pour le transmultiplexeur complet. 
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Additions reelles 


Multiplications 
reelles 


Realisation transverse 
Realisation en echelles 
Realisation en treillis 
(normalise) 


2m + 31og 2 2M 
4m + 4 + 3 log 2 2M 
4m + 3 Iog 2 + 2M 


2m + 4 + 2 log 2 2M 
4m + 6 +2 log 2 2M 
4m + 4 + 2 log 2 2M 



TABLEAU 2 - Nombre d'operations reelles par sous-bande et par echantillon 
pour le modulateur (ou le demodulateur). 



Le gain apporte par rapport a une realisation directe du schema de la 
figure 1 est done net, puisque celle-ci necessiterait 2mM - 1 additions complexes 
et 2mM + 1 multiplications complexes par sous-bande et par echantillon, au 
modulateur comme au demodulateur. 

En termes de cases memoire, il faut stocker 4M valeurs complexes pour 
realiser la pre-modulation ainsi que les coefficients des differentes structures. 
Lorsqu'on a les memes filtres a remission et a la reception, on obtient la 
premiere colonne du tableau 3. Par ailleurs, il faut dans tous les cas stocker 
4(m + i)M valeurs complexes dans un "buffer" pour le filtrage polyphase au 
modulateur comme au demodulateur. 





ROM 


RAM 


Realisation transverse 


2(m+ l)M+2 


4(m + 3) M 


Realisation en echelles 


(2m + l)|^^-i + 2M + 2j 


4(m + 3) M 


Filtre transverse 
symetrique 


(m + 2) Af+ 2 


4(m + 3) M 


Realisation en treillis 
1 (normalise) 


m|^yJ+2M + 2 


4(m + 3) M 



TABLEAU 3 - Cases memoire reelles pour le modulateur (ou le 
demodulateur) complet. 
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Les differentes techniques proposees se caracterisent notamment par le 
fait que, pour un systeme modulateur-demodulateur mis "dos-a-dos", leurs EES 
et IEC sont exactement nulles. En pratique, du fait de 1'imprecision du calcul 
numerique, elles sont generalement de l'ordre de 10* 14 . 
5 Dans le cas des modes 2 et 4, cette caracteristique de reconstruction 

parfaite est assuree structurellement, c'est-a-dire qu'elle est maintenue apres 
quantification des coefficients echelles pour le BFDM/OM ou treillis pour 
TOFDM/OM. 

Deux criteres peuvent etre pris en compte pour le design des filtres 
10 prototypes : la localisation et la selectivity On peut egalement tenir compte 

d'autres aspects, tels que les distorsions de canal representatives de differents 
canaux de transmission, par exemple du type radio-mobile. 

A titre d'exemples purement indicatifs, les tableaux 4 et 5 de l'annexe D 
donnent des realisations particulieres de Invention, dont les resultats sont 
15 illustres par les figures 1 1 A, 1 IB, 12A et 12B. 

Les figures 11A et 11B presentent respectivement la reponse temporelle 
et la reponse frequentielle pour un prototype biorthogonal avec M - 4, N = 32, 
a= 8, § = 0,9799 (localisation), 0,9851 (localisation modifiee, selon le 
critere de Doroslovacki). Elles correspondent a la premiere colonne du tableau 4 
20 (coefficients transverses) et au tableau 5 (coefficients des echelles). 

Les figures 12A et 11B presentent respectivement la reponse temporelle 
et la reponse frequentielle pour un prototype biorthogonal avec M = 4, N = 32, 
a= 2, £ = 0,9634 (localisation), ^ = 0,9776 (localisation modifiee, selon la 
mesure de Doroslovacki). Elles correspondent a la seconde colonne du tableau 
25 4. 
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Modulation mult iporte use de type 
BFDM/OQAM 

Dans cette annexe, en guise d'introduction aux modulations BFDM/ OQ AM, 
nous rappelons quelques definitions essentielles sur la biorthogonalite ([16], 
[17], [18]). 

Soit E un espace vectoriel sur un corps K, les definitions et propnetes que 
nous allons utiliser pour generer une modulation BFDM/OQAM peuvent se 
resumer ainsi : 

Definition A.l Soient (xi) ie i et (x;) l€ i deux families de vecteurs de E. (x,), € i, 
(it)*€i sont biorthogonales 

si et seulement si: € I 2 , (xi 7 Xj) = S{j 

Definition A.2 Soient (xi) ie i et (x, ) l€ i deux families de vecteurs de E. (x t -) t€ i, 
(xi)iei forment un couple de bases biorthogonales de E si et seulement si: 

• (xijiei et (xi)iei forment deux oases de E 

• (#;)iei e * (zi)iei soni deux families biorthogonales 

Propriete A.l Soit ((x t )* € i> (*i)iei) un couple de bases biorthogonales de E, 
alors Vx G E: 

• x = 52 ie i(xi,x)xi = J2iei{*i> x ) x i 

• 5t x = Y^iei a » x «V a ^ rs a * = x ) 

• si x = ^2 ie i &iXi, alors 6c{ = (x;, x) 

• IN 2 = E.-€i(^*>'< 5 ^ x > 

Un signal complexe module en frequence sur 2M sous-porteuses peut s'ecrire 

+00 2M-1 

S (t)= £ - nro)e 2j ' ir(/o+m " o) V' v '"- (1) 

71= — CO 771=0 

avec: 

_ Gm,n € R; 

-Aun filtre prototype reel, de largeur de bande v 0 et de support fini: 
h(t) e'[-T u T 2 } avec T x et T 2 des reels; 
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- fo = 0; 

Pour obtenir une modulation biorthogonale on cherche a ecrire s(t) a l'aide 
dun couple (Xm^Xm.n) de bases biorthogonales : 

+00 2M-1 

= 53 S fl m>nXm,n(0 ( 2 ) 

avec: 

a m ,n = 5? {y^3(*)X»,n(0*J = (^Xm,n> (3) 

La derivation des expressions des bases discretes associees est presentee en 
annexe 2. i 

Apres translation de 7\ et discretisation a la periode T e = j$ = 2M^' ^ est 
egalement possible, cf. annexe 2, de definir un couple de bases biorthogonales 
discretes (xm,nW, Xm,n[t]) tel que 

H-oo 2M-1 

S W = 53 53 a m f nXm,n[*] W 
n=s— oo m=:0 

avec: 

Xm,n[k] = (-l) mn v^ p[* - n M]e J '^ m ( fc - nM -^ M) (6) 
*»»!»[*] = (-lr^gfAr-nMy^*-*^-^ (7) 
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Annexe B 
Le transmultiplexeur BFDM/OQAM 

B.l Cas general biorthogonal 

On pose N la longueur du filtre prototype p[k], de telle sorte que: 

2T = T 1 +T 2 = (N - l)T e (8) 
et Ti = 2AT, T 2 = 2(1 - A)T avec A € [0, 1]. Mors: 

+oo 2M-1 

s[k] = £ a™A-V mn Plk-nM} 

n=— oo m=0 

x e J'l& m (*-'» w )e i ( Vm - n -i& mA( ^" 1) ) (9) 

+co 2JW-1 
ti— — oo m=0 

xei^-C^-^JeiC^+^^C^-^^" 1 ))) (10) 

avec Z> un parametre fixe arbitrairement et qui, comme on le verra, permet de 
gerer le retard de reconstruction. Au vu de l'equation (10), on pose maintenant: 

x m (n) = (-l) m V' ?n ^ (12) 
/ m (Jfc) = V^pfifc]^^*-^ 1 ) (13) 

de sorte que: 

+oo 2M-1 

s[fc] = £ J] x m (n)/ m (*-nM) ( 14 ) 

tj.=— oo m=0 

Pax ailleurs, la base duale de demodulation s'ecrit: 

X m ,n[k) = (-ir n v^ «[* - nJfle'aKw-Wf (15) 
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a mn = &{(-l) mn v/22>[%^- nM ] 

k ^ 

xe ->m< k - n ' *- J We'- J "*"> ( 16 ) 

= sJ(-l) mB \/5]£s[£> + nA/- %[!>-*] 
' 1 



(17) 



(18) 

(19) 
(20) 

(21) 



ce qui nous amene a poser: 

h m (k) = V2q[D- fc^*"**"***) 

et 

D _ a M _ p avec a et 0 entiers etO</?<M-l 
de sorte que: 

a m n-a = » |(-l) m ( n - a) e- i ? (n ~ Q) 52 s ^ nM ~ k ~ fl hm ( k ^ 

Le facteur (-1)™ apparait aussi bien au modulateur qu au modulateur si bien 
qu'on peut ■ le supper sans rien changer, et on aboutit alors au schema du 
transmultiplexeur de la figure 1. 

B.2 Cas particulier orthogonal 

Dans le cas orthogonal, on a D = N - 1 et q[k] = p[k], d'ou: 

Mfc) _ V2 P [iV-l-fcy^-^e-?- (23) 

Dans le cas ou <*[*} = p[D - k\ on a alors rffc] = P [N - 1 - P^£«? 
svmetrique. Mais contrairement a ce qu'on peut souvent lore d ^X!ument 
cite ou explicite ([4], [6], [7], [9]), la symetne du prototype n est absolument 
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pas necessaire. On pourra pour s'en convaincre prendre Tun des prototype sui- 
vants et verifier numeriquement (une verification directe est assez fastidieuse) 
qiril permet la reconstruction parfaite pour M = 4 dans le cas orthogonal (on 
a alors A r - 1 = D = 7, a = 2 et 0 = 1): 

P(z) = i (1 + z- x + *~ 2 + z- 3 - z- 4 - z- A - z- 5 - z- 6 - z- 7 ) (24) 
p(-) = l(l + 2 -i-z-* + z- 3 + z- 4 + z- 4 + z- 5 -z-* + z- 7 ) (25) 

On peut meme verifier que tout prototype P{z) = ]CL=o P( n ) z ~ n verifiant 



(26) a (31) assure aussi la reconstruction parfaite pour M = 4 dans le cas 
orthogonal: 

Co = ±1 , d = ±l (26) 

e 0 = ±l , £i = ±l (27) 

W0)|<^ , b(i)l<x (28) 

p{ 2)=e iy j\- P {iy , P(3)=eo^/|-P(0) 2 (29) 

p(4)=co £ o-y/i-p(0) 2 , p(5)=c l£lV /i-p(l) 2 (30) 

p(6) = Cl £xKl) , P(7) = OBffuKO) ( 31 ) 
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Annexe C 
La condition de biorthogonalite 

C.l Approche polyphase 

La realisation cTun schema de modulation et demodulation selon la figure 
1 serait extremement couteux en termes de complex^ <*°^£J%f°T 
position du prototype P(z) en fonction de ses composantes polyphases Gj(z), 

telle que 



2M-1 



(32) 



z=o 



permet d'exprimer ensuite les bancs d'analyse et de synthese sous la forme 



1=0 
2M-1 



41V1 — " X 

1=0 



(34) 



On en deduit ainsi l'expression des matrices polyphases R(z) et E(z) des 
bancs de filtres du modulateur et demodulateur : 



E(z 7 ) = Wi 



i*(z 2 ) = J 




(35) 



G 2M -i(z 2 ) 
G 0 (z*) 0 

0 G 2M -i(z 2 ) 
G 2 m-i{* 2 ) 



Wl 



(36) 



<?o(z 2 ) 0 
ou J, et W 2 sont definis ci-dessous: 

0 1 



J = 




(37) 
(38) 
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De maniere a isoler la fonction transmultiplexeur nous introduisons lajiota- 
tion qui suit, ou par commodite d'ecriture, on ne prendra plus x m (n) = j n a m , n , 
mais x m (n) = a m ,„, de maniere a avoir: 



/ Xm(n) = a m , n ^ | Xmi-J 2 

{ R(x' m (n)) = x m (n) = a m , n \ X' m {-jz 



-jz) i — > J R a m ,„ 



(39) 

x m (n) represente les symboles reels a emettre et x' m (n) les symboles complexes 
regus avant extraction de la partie reelle. La figure 2 donne une vue globale de 
la chaine. 

Les matrices polyphases B(z 2 ) et R{z 2 ) permettent d'obtemr une repre- 
sentation sous forme polyphase du transmultiplexeur (figure 3). H reste alors a 
prendre la partie reelle des echantillons de sortie x' m (n - a) pour reconstituer 
l'entree avec un retard de a echantillons. 



C.2 Relation d'entree-sortie 

Nous notons X{z) le vecteur representant, dans le domaine transforme 
en z, les donnees emises. A la reception, apres demodulation, nous notons 
15 par x'(-jz) le vecteur des transformees en z associe aux donnees regues. 

L'extraction de la partie reelle fournit ensuit le vecteur X(z). Notre but est 
alors : 

- de determiner la relation entree-sortie, c'est-a-dire la relation entre X(z) 
etX(z); 

- de determiner les conditions sur les composantes polyphases G t (z) de 
20 permettant de garantir l'egalite X(z) = X(z) ; 

- d'en deduire le retard de construction a. 

Les 3 principaux elements de ce schema permettant de determiner la re- 
lation entree-sortie sont les 2 matrices polyphases E(z) et R(z) ainsi que la 
matrice de transfert A 0 {z), liee aux expanseurs, delais et decimateurs. Pour 
determiner cette derniere on peut se baser sur le cas elementaire represente a 
25 la figure 4 pour lequel la fonction de transfert est donnee par 

JO si A' n'est pas multiple de M, ^ 4( ^ 

V ^ = \ z-% U{z) si K est multiple de M 
De la figure 3 il vient ensuite : 

z -* X \-jz) = E{z 2 )A 0 {z)R{z 2 )X{-jz) 
(jz)- a x'(z) = E(-z 2 )Ap(jz)R(-z 2 ))X(z) 

X\z) = j a z a W l G(jz)WlX(z) (41) 
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oil la matrice G[z) est definie ci-dessous : 



G{z) = 
V 0 



0 \ / o 



G2M-l(^ 2 ) 

0 



On a alors: 

X(z) = Q(*)X(*) 

avec: 

Apres calcul, on obtient: 



(42) 
(43) 
(44) 



Q{z) = 



( Q 0 (z) 0 Qi(z) ••• Q*f-iW 
0 Qo(z) 0 <?i(*) 
0 Qo(z) '•• 



<?i(*) 



0 

V 0 Qm-i{z) ■■■ Qi(*) 0 



0 \ 
Qm-i( z ) 

0 

Qo(*) J 



avec: 



Q 4 (z) = 2(-l)< Y, U < (~ z ) COS V " 2 J J 

J=0 



(45) 



(46) 



La signification de d est precisee plus loin. L'expression exacte de M~z*) 
depend d'ailleurs de ce parametre d (entier positif ou mil), on pent alors de- 
montrer qu'on a reconstruction parfaite si et seulement si: 



si 0 < d < M - 1: 

- si 0 < I < d: 

z~ s 

Gt(z) G d -i{z) + z~ l G M +i{z) G M +d-i(z) = 2^ 

- si d + 1 < I < M - 1: 

G,(z) GW+«i-t(*) + G M +i(z) G M +d-i(z) = ■ 2 jjf 



(47) 



(48) 
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- si M < d < 2M - 1: 

- si 0 < / < d - M: 

G { (z) G d .i{z) + G M+ i{z) G d - M -i{z) = ( 49 ) 

- si d + 1 - M < I < M - 1: 

G t (z) G d .i{z) + z~ l G M+l {z) G M +d-i(z) = (50) 



dets les entiers definis par D = 2sM + <*, 5 > 0 et 0 < d < 2M - 1. Le 
retard de reconstruction a est lie au parametre s par les relations: 



avec 



'2s si d = 0 

a= i 23 + 1 sid e{l,...,M} (51) 
2(5 + 1) sid€{M + l,...,2M-l} 

De ce resultat on peut deduire le cas particulier orthogonal pour lequel D = 
N - 1, avec N la longueur du filtre prototype paraunitaire, c'est-a-dire sy- 
metrique ici (on dit que P{z) est paraunitaire si P(z) = z~^ N -^P(z) avec 
15 P{z) =■ P'(z~ 1 )). On peut en effet verifier que: 

G d -i{z) = z-'Gi{z) . si Q < / < d (52) 
G 2M+d -i = z- (3 - l) G,(z) si d + 1 < / < 2M - 1 (53) 

Ainsi, dans le cas particulier orthogonal, on a reconstruction parfaite, avec un 
retard a = N ~m +P ■> S1 et seulement si: 

20 Gi(z) G t (z) + G l+M (z) G M+l (z) = ^ 0 < / < M - 1 (54) 



25 



30 
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Annexe D 

E Coefficients des nitres prototypes obtenus 
par optimisation 



o 
i 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 



exemple figures 11 A, 11B 



-4.460868105953324e-05 
-8.827693704913472e-05 
1. 816721975 145588e-04 
2.302861759368 1 1 le-04 
-1.172447599636272e-03 
-1.827055790732281e-03 
-3.760044826875730e-03 
-6.052599354959620e-03 
-6.541334278009250e-03 
-2.320957844665448e-03 
1.063973252261601e-02 
4.151536856291601e-02 
1 .043838059706333e-01 
2.005189128209921e-01 
2.913131449163113e-01 
3.352627462532674e-01 
3.351172696026857e-01 
2.909415993397522e-01 
2.000454638421703e-01 
1.039959799288574e-01 
4.124129545474275e-02 
1.040407270162191e-02 
-2.440939909195805e-03 
-6.649467765409867e-03 
-5.878063999471562e-03 
-2.983359331348727e-03 
-2.475652683945518e-03 
-8.249232623623326e-04 
-8.891453128240245e-05 
-4.464223698699074e-04 
3.70450426982971 le-04 
1.247820119405080e-05 



exemple figures 12A, 12B 



5.014949968230972e-02 
1 .4555838 164890 19e-0 1 
2.500737066757044e-01 
3.228805982747062e-01 
3.661515202615520e-01 
3.515099229023638e-01 
2.548572622632939e-01 
1.351315191913547e-01 
5.168757567424829e-02 
1 . 125177964848224e-02 
-4.802381010162210e-03 
-1 . 106221296463278e-02 
-8.872655589434630e-03 
-3.753426678003194e-03 
-1.654867643757010e-03 
-1.383187971152199e-03 
-3.813932123570836e-04 
5.624569918905843e-06 
-2.475635347949224e-06 
4.952234305253537e-05 
2.048787314180216e-05 
1.004915628115845e-07 
4.423431446835775e-08 
2.649149072234273e-06 
1.785190585201287e-08 
6.863128678632993e-12 
-3.021006714034441e-12 
-2.308376067571482e-09 
6.598349790230041e-l2 
-5.789492272166598e-16 
-2.548414264695020e-16 
8.532126971482393e-13 



TAB. 4 - Prototypes biorthogonaux avec M = 4 et N = 32 (coefficients 

transverses). 
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1=0 


1=1 


/o,o 
/0.1 

/o,2 

eg 

68 
*S 
*6 

4 


2.S20843813510179e-01 
-5.673498928902276e-01 

3.721645241266496e-01 
-1.1 67480680205269e-01 

1.064716331427927e-01 
-1.596131503239696e-01 
-2.626782049187459e+01 
- 1 .61444704546251 le-04 
2. 628294699 12271 7e+01 


1.299280891559943e-01 
-3.310904763283146e-01 
-3.738170157940610e-02 
-2.847704852297620e-02 
-1.181506919769764e-01 

9.127789670781582e-02 
-2.054722686392923e+01 
-3.551599154933042e-04 
2.069686434312222e+01 



TAB. 5 - Prototype biorthogonal (coefficients des echelles) avec M = 4, 
10 N = 32 et a = 8 (cf. fig. 11A, 11B). 
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REVENDICATIONS 

1. Procede de transmission d'un signal mulliporteuse biorthogonai BFDM/OM, 
caracterise en ce qu'il met en oeuvre une structure de transmultiplexeur assurant 

5 : 

une etape de modulation, a l'aide d'un banc de filtres de synthese 

(11) , presentant 2M branches paralleles, M > 2, alimentees 
chacune par des donnees source, et comprenant chacune un 
expanseur d'ordre M et des moyens de filtrage ; 

10 - une 6tape de demodulation, a l'aide d'un banc de filtres d' analyse 

(12) , presentant 2M branches paralleles, comprenant chacune un 
decimateur d'ordre M et des moyens de filtrage, et delivrant des 
donnees re§ues representatives desdites donnees source, 

lesdits moyens de filtrage etant deduits d'une fonction de modulation prototype 
15 predetermine. 

2. ProcedS de transmission selon la revendication 1, caracterise en ce que lesdits 
moyens de filtrage dudit banc de filtres de synthese et/ou dudit banc de filtres 
d' analyse sont respectivement regroupes sous la forme d'une matrice polyphase. 

3. Proc6de de transmission selon la revendication 2, caracterise en ce qu'au 
20 moins une desdites matrices polyphases comprend une transformee de Fourier 

inverse(51, 61) a 2M entrees et 2M sorties. 

4. Proced6 de modulation d'un signal transmis selon le procedd de l'une 
quelconque des revendications 1 k 3, caracterise en ce qu'il met en ceuvre une 
transformee de Fourier inverse (51) alimentee par 2M donnees source ay ant 

25 chacune subie un decalage de phase prdd&erminee, et alimentant 2M modules 

de filtrage, suivis chacun d'un expanseur d'ordre M, dont les sorties sont 
regroupees puis transmises. 

5. Procede de modulation selon la revendication 4, caracterise en ce qu'il 
delivre des donnees s[k] telles que : 

30 x° m (n) = a m y* n 

2M-1 1* D-M 2n_ 



x l t (n) = j2^x° k (n)e '™ 2 e*'< 
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D-M\ 



2 M -JllFFTi x° (n), jcJ„_, (n)e 



z [I] 



m-i 



xf (n) = £ P(/ + 2*Af ) (w - 2*) 




et 



ou 



avec 



D = aM - P, 

a entier representant le retard de reconstruction ; 
P entier compris entre 0 et M-l ; 
|_.J est la fonction "partie entiere". 



6. Procede de demodulation d'un signal transmis selon le procede de Tune 
quelconque des revendications 1 & 3 caracteris£ en ce qu'il met en ceuvre une 
transformee de Fourier inverse (61) alimentee par 2M branches, elles-memes 
alimentees par ledit signal transmis, et comprenant chacune un d€cimateur 
d'ordre M suivi d'un module de filtrage, et alimentant 2 M multiplieurs de 
ddcalage de phase, ddlivrant une estimation des donnees source. 
0. ProcedS de demodulation selon la revendication 6, caracterise en ce qu'il 
ddlivre des donnees a^^ telles que : 



(n-a) = s[nM -p-l] 



j?J (n - a) = ]>>(/ + 2kM)2]{n -a -2k) 





= 2Mj2e~ ,2M 2 



IFFT(jc;'(, 



n 




avec : D = 2.S.M + d, 
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oil : s est un entier ; 

d est compris entre 0 et 2M-1. 

8. Procede de modulation selon Tune quelconque des revendications 4 et 5 ou 
de demodulation selon Tune quelconque des revendications 6 et 7, caracterise en 
ce que lesdits modules de filtrage sont realises sous Tune des formes 
appartenant au groupe comprenant : 

les filtres a structure transverse ; 
les filtres a structure en echelle ; et 
les filtres a structure en treillis. 

9. Procede selon Tune quelconque des revendications 1 h. 8, caracterise en ce 
que ledit signal multiporteuse biorthogonal est un signal OFDM/OM. 

10. Dispositif d* emission et/ou de reception d'un signal BFDM/OM, 
mettant en ceuvre le procede de Tune quelconque des revendications 1 a 9. 
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